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A computational strategy and a collection of codes are presented for studying multiparticle final states in hard 
hadronic collisions. 



1. Introduction 

Multijet final states are characteristic of a large 
class of important phenomena present in high- 
energy hadronic collisions. QCD interactions 
generate multijet final states via radiative pro- 
cesses at high orders of perturbation theory (pure 
QCD processes). Heavy particles in the Stan- 
dard Model (SM), such as W and Z bosons or 
the top quark, t, decay to multiquark configura- 
tions (eventually leading to jets) via electroweak 
(EW) interactions (mixed EW-QCD processes). 

In addition to the above SM sources, particles 
possibly present in theories beyond the SM are ex- 
pected to decay to multiparton final states, and 
therefore to lead to multijets. Typical examples 
are the cascade decays to quarks and gluons of su- 
persymmetric strongly interacting particles, such 
as squarks and gluinos (SUSY-QCD processes). 

In addition to fully hadronic multijet final 
states, a special interest exists in final states 
where the jets are accompanied by gauge bosons. 
For example 
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Z + jets — ► vv + jets 



(1) 



is an important background to SUSY searches. 
Likewise, 
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W + 4 jets — > t v t + 4 jets 



(2) 



provides the leading source of backgrounds to 
the identification and study of top quark pairs 
in hadronic collisions. 
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Several parton-level Monte Carlo (MC) event 
generators exist in the literature (such as VECBOS 
[0], MADGRAPH |§, CompHEP g, GRACE [Q, 
PAPAGENO ||) and studies can be performed by 
assuming that hard partons can be identified with 
jets, and that the jets' momenta are equal to those 
of the parent partons. This simplification gives 
rough estimates, but cannot be used in the con- 
text of realistic detector simulations, for which 
a representation of the full structure of the final 
state (in terms of hadrons) is required. 

This full description can be obtained by merg- 
ing the partonic final states with shower MC pro- 
grams where partons are perturbatively evolved 
through emission of gluons, and subsequently 
hadronized. 

However, this merging is not always possible, 
since common parton-level MC's sum and aver- 
age over flavours and colours, and do not usually 
provide sufficient information on the flavour and 
colour content of the events. 

In this contribution we review a strategy for 
the construction of event generators for multi- 
jet final states, based on the exact leading-order 
evaluation of the matrix elements for assigned 
flavour and colour configurations, and the sub- 
sequent shower development and transition into 
a fully hadronized final state. 

We also present ALPGEN || , a collection of codes 
realized within the framework of the presented 
approach. 

As a case study, we then focus on the produc- 
tion of WQQ + n jet final states [0 (with Q being 
a massive quark, and n < 4) and present results 
for several production rates and distributions of 
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interest at the LHC and at the TEVATRON. 



3. The ALPGEN project 



2. The computational strategy 

There are two main sources of problems: the 
matrix-element evaluation and the reconstruction 
of the colour flow. To illustrate our solutions we 
take, for simplicity, the case of multigluon pro- 
cesses, as the extensions to cases with quarks and 
EW particles follow the same pattern. 

The process gg — > 10 g has 5 x 10 9 contributing 
Feynman diagrams. To efficiently carry out the 
computation, an algorithm is therefore needed 
which does not make explicit use of Feynman di- 
agrams. Fortunately, such an algorithm exists 
||||. The key ingredient is the use of a recur- 
sive procedure to determine the matrix element 
directly from the effective action. 

To better illustrate the problem with the color 
flow, we take again the case of multigluon produc- 
tion. The scattering amplitude for n gluons with 
momenta pf , helicities e^ 1 and colours <2j (with 
i = I, ... ,71), can be written as 

M({ Pi }, M, {a t }) = £ MA - A a - . . . A a - ) 

P(2,3,...,ra) 

A({ Pil },{e il };...{p i J,{e in }), (3) 

where the sum extends over all permutations Pi 
of indices (2,3,..., n), and the functions A{{Pi}) 
(known as dual or colour- ordered amplitudes) are 
gauge invariant, cyclically-symmetric functions of 
the gluons' momenta and helicities. 

Each dual amplitude A({Pi}) corresponds to 
colour flows from one gluon to the next, accord- 
ing to the ordering specified by the permutation 
of indices. Furthermore, each different colour or- 
dering corresponds to a different shower evolution 
of the event, once the hadronization programs are 
turned on. Therefore, any realistic event gener- 
ator should keep track, event by event, of the 
colour flow (and flavour content) to perform a re- 
liable simulation of the subsequent evolution into 
fully hadronized final states. 

More details on our algorithm for the colour 
reconstruction can be found in ref. 10 . 



ALPGEN m is a collection of codes for the gen- 
eration of multi-parton processes in hadronic col- 
lisions implementing the computational strategy 
described in the previous section. 

The processes currently available in the pack- 
age are: 

^ -» WQQ + n jets (n < 4) 

-> Z/j*QQ + n jets (n < 4) 

-> W + n jets (n < 6) 

-> Z/j* + n jets (n < 6) 

-> n z Z + n w W + rih H + n jets 
n w + rih + n < 8, n < 3) 

-> QQ + n jets (n < 6) 

-> QQQ'Q' + n jets (n < 4) 

-> HQQ + n jets (n < 4) 
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Q and Q' being massive quarks. Anywhere, ex- 
cept for processes 5), W and Z/^f* are consid- 
ered off-shell, meaning that the corresponding 
produced final state is a pair of leptons. 

The documentation and the codes can be found 



at http : //mlm . home . cern . ch/mlm/ alpgen/ 



4. Results 

We show a few illustrative examples obtained 

with ALPGEN for the process p'p 1 —* Wbb + n jets. 
At the LHC one gets the partonic rates of Table |l|. 
The numbers in the first column refer to the sub- 
processes given in ref. . The qualitatively new 
processes appearing when Nj increases are large 
and are responsible for the growth of the partonic 
rates when going from Nj = 2 to Nj = 3. 

In Figure [l] we show a comparison between the 
pr distributions of jets at the TEVATRON before 
and after the perturbative phase of the shower 
evolution performed by HERWIG pl[ . A good 
matching of the jet spectra is achieved when us- 
ing generation (parton level) cuts loser than jet- 
defining cuts. Finally, we address the issue of the 
ability of the shower MC to correctly predict the 
rate for hard radiation leading to extra final-state 
jets. In Figure we compare the jet rates eval- 
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Table 1 



Partonic rates in pb at the LHC, as a function of the total number of jets Nj (including b and b). 



Process 


Nj = 2 


Nj = 3 


Nj = 4 


Nj = 5 


Nj = 6 


1 


2.60(1) 


0.63(1) 


0.144(3) 


0.036(2) 


0.008(1) 


2 




2.97(1) 


2.11 (2) 


1.08(2) 


0.47(2) 


3+4 






0.288(1) 


0.24(1) 


0.13(2) 


5 








0.030(1) 


0.031(4) 


6 










0.0010(3) 


Total 


2.60(1) 


3.60(1) 


2.54(2) 


1.38(2) 


0.64(3) 


The jet defining 


cuts are p l T > 20 GeV, \rn\ 


< 2.5 and ARij 


> 0.4. 
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Figure 1. Inclusive pt distributions of jets at the 
parton level, with separation cut AR > 0.7 (solid 
curves), and of fully-showered Rj et — 0.7 jets 
(dashed curves). These last are obtained starting 
from a full sample of AR > 0.4 partonic events. 



uated with exact matrix elements with those ob- 
tained from hard radiation during the shower evo- 
lution of lower-order parton-level processes. We 
find a good agreement in the range pr ^ 45 GeV, 
for the radiation of one extra jet by HERWIG, while 
not enough hard radiation is emitted by HERWIG 
to correctly predict the emission of 2 extra jets. 

5. Summary 

The formalism we have outlined allows us to 
compute several complex multiparton processes. 

ALPGEN is a collection of codes realized within 
the framework of the presented approach. It can 



E I4 m f + 4 jets 
Soiid: W+4j ME (PR0C=l+3) (22 fb) 
Dashes: W+4j ME + Herwig (13 fb) 
X: T+3j ME + Herwig (15 fb) 

W+2j ME (PR0C=1) + Herwig (6 fb) 
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Figure 2. p T distributions of jets at the TEVATRON: 
Matrix Element versus Parton Shower. 



be used to perform realistic studies of multiparti- 
cle processes at TEVATRON and LHC. 

ALPGEN is ready-to-use and available at 
http//mlm.home . cern. ch/mlm/alpgen/. 
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